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摘 要 : 基于 CIMMS NDVI 3g v1.0 数据 集 和 日 值 气象 数据 ,结合 极端 气温 指数 , 辅 以 极点 对 称 模 态 分 解 .趋势 分 


BFF Mann-Kendall 趋势 检验 .相关 分 析 等 方法 ,探讨 中 国 北方 生长 季 植 被 覆盖 及 极端 气温 的 变化 特征 ,研究 植被 覆盖 
对 气温 极 值 的 响应 状况 。 结 果 表 明 :GD 1982—2015 年 中 国 北方 生长 季 NDVI 以 0.002 + (10a) ”的 速率 上 升 (P< 
0. 05) ,ESMD( 极 点 对 称 模 态 分 解 方法 ) 显示 生长 季 NDVI 波动 上 升 ; 针 叶 林 , 灌 从 .荒漠 植被 .草地 以 及 栽培 植被 时 
增长 趋势 ,栽培 植被 增 速 最 快 , 针 冰 混 交 林 、 沙 叶 阔 叶 林 和 高 山 植被 星 不 显著 减少 趋势 。G@ 空间 上 ,NDVI 显著 增加 
区 域 超过 全 区 的 33% ,主要 分 布 在 天 山 .塔里木 盆地 北部 .祁连山 . 院 南 山区 .黄土 高 原 .河套 平原 . 昌 梁 山 和 太行 


山 \ 大 别 山 以 及 辽西 丘陵 地 区 ;显著 下 降 区 域 仅 占 12% ,主要 分 布 在 大 兴安 岭 .小 兴安 岭 和 长 白山 区 。@) 极端 气温 


指数 中 , 除 TNmean( 日 最 低 气温 平均 值 ) 和 TN 日 最 低 气 温 极 低 值 ) 呈 上 升 趋势 外 ,其 余 冷 极 值 指数 均 呈 下 降 趋 


势 ; 所 有 暖 极 值 指数 均 呈 上 升 趋势 ;其 他 指数 中 ,DTR( 气 


温 日 较 差 ) 呈 减 小 趋势 ,GSL( 生 长 季 日 数 ) 呈 增 加 趋势 。 


@ 中 国 北方 NDVI 与 极端 气温 指数 的 相关 性 表明 , 冷 极 值 指数 中 NDVI 与 FDO (霜冻 日 数 ) 、TN10p( 冷 夜 日 数 )、 
TX10p( 冷 昼 日 数 ) 呈 显著 负 相 关 (P<0.05) ,与 TNmean 呈 显 著 正 相关 (P<0.01) ;NDVI 与 所 有 了 暧 极 值 指数 呈 正 相 
关 , 与 TR20( 热 夜 日 数 ) .TXmean( 日 最 高 气温 平均 值 ) .TX90p( 暧 慎 日 数 ) 以 及 TN90p( 暖 夜 日 数 ) 存 在 显著 相关 性 
(P<0.05);NDVI 与 GSL 旦 显著 正 相 关 (P<0.05)。@ 天山、 塔里木 使 地 北 缘 、 祁 连 山区 河套 平原 、 黄 土 高 原 、 太 
行 山 和 扣 梁 山区 等 NDVI 显著 增加 区 域 对 极端 气温 指数 的 响应 强烈 。NDVI 显著 增加 区 主要 对 FD0、TNmean、 
TN90p、GSL 等 指数 响应 较 强 。NDVI 显著 减少 区 域 对 指数 的 响应 各 异 , 主要 与 SU25 (夏季 日 数 ) 呈 显著 负 相 关 


(P<0.05)。 


关键 词 : 生长 季 ; GIMMS NDVI 3g v1. 0 数据 集 ; 时 空 变化 ; 极端 气温 指数 ; 响应 ; 中 国 北方 


气候 变化 正在 影响 全 球 生 态 系统 ,植被 作为 全 
球 生态 系统 的 重要 组 成 部 分 ,对 气候 变化 尤为 敏感 。 
近年 来 ,确定 气候 变化 如 何 影 响 植被 变化 受到 了 全 
BREA WI IZ ET. AJE RPK TS LA Te 
漠 生 态 系统 等 均 是 对 气候 变化 敏感 的 生态 区 
a9) 。 了 解 植被 动态 变化 并 确定 其 对 气候 变化 的 
响应 情况 ,对 确定 并 保持 生态 系统 的 质量 和 状态 至 
KEAT, 

气候 变化 分 为 波动 变化 .趋势 变化 以 及 极端 事 
FZX., IPCC 第 五 次 气候 变化 评估 报告 指 
出 ,气候 变 暖 不 仅 包括 平均 状态 的 气候 变化 ,还 包括 
极端 气候 变化 " 。 过 去 一 百年 中 国 地 表 年 均 气温 
呈 显 车 上 升 趋势 ,北方 明显 变 暖 ,极端 气温 事件 增 
ZO). XB HEE 指出 ,中 国 北方 气温 上 升幅 度 
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高 于 全 国平 均 气 温 ,干燥 和 变 暖 是 区 内 气候 变化 的 
突出 趋势 。 世 界 气象 组 织 (WMO) 气候 学 委员 会 
(CCL) ,全 球 气候 研究 计划 (WCRP) 气候 变化 以 及 
可 预测 性 计划 (CLIVAR) 气候 变化 检测 监测 和 指 
标 专家 组 (ETCCDMI) 制定 的 极端 气候 指数 ,被 用 于 
极端 气候 的 分 析 。 众 多 学 者 基于 极端 气候 指数 对 桔 
阳 湖 流域 ,中国 天 山南 北 坡 \ 呼 伦 贝 尔 草 原 、 友 岭 山 
地 陕西 段 等 地 区 的 极端 气候 ,特别 是 气温 极 值 变化 
情景 进行 了 研究 “" 。 除 了 关注 极端 气候 事件 本 
身 外 ,近年 来 也 有 学 者 分 别 对 新 疆 、 陕 甘 末 地 区 、 西 
北 干 旱 区 以 及 治 海地 区 的 极端 气候 事件 对 植被 覆盖 
变化 的 影响 进行 了 深入 研究 。 这 些 研究 主 
要 使 用 线性 趋势 分 析 法 对 NDVI 及 气候 因子 的 时 间 
变化 趋势 的 分 析 ,未 能 较 好 地 反映 要 素 在 各 阶段 的 
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变化 ,在 变化 程度 上 部 分 夸大 或 缩小 了 要 素 的 上 升 
HBO”, 2013 年 极点 对 称 模 态 分 解 方法 (ESMD ) 
的 提出 ,为 分 析 不 同 要 素 时 间 序 列 的 非 线 性 趋势 
变化 特征 提供 了 新 的 研究 思路 。 

气候 变 暖 背景 下 ,Sun 4) 基于 GIMMS NDVI 
数据 对 1982 一 2006 年 中 国 北 方 的 植被 动态 进行 评 
估 后 ,表明 中 国 北方 植被 覆盖 率 呈 小 幅 上 升 趋势 , 且 
草地 和 农田 的 NDVI 增加 较 多 。 中 国 北方 不 同 区 域 
的 植被 覆盖 也 发 生 了 明显 变化 下 。 这 些 研究 还 
显示 ,气温 对 植被 变化 的 影响 强 于 降水 。 基 于 此 背 
景 ,本 研究 利用 最 新 发 布 的 CIMMS NDVI 3g v1.0 数 
据 和 日 值 气象 数据 ,选取 极端 气温 指数 ,运用 趋势 分 
析 、 极 点 对 称 模 态 分 解 、Mann-Kendall 趋势 检验 、 相 
关 分 析 等 方法 ,对 中 国 北方 植被 覆盖 时 空 变化 特征 、 
极端 气温 变化 特征 以 及 植被 覆盖 与 气温 极 值 之 间 的 
关系 进行 分 析 , 这 对 了 解 近 34 a 中 国 北方 植被 覆盖 
变化 和 气温 极 值 对 植物 生长 的 影响 、 推 进 中 国 北方 
生态 工程 建设 和 保障 生态 安全 方面 有 重要 意义 。 


1 数据 与 方法 


1.1 研究 区 概况 

中 国 北 方 (31?23' ~ 53°31'N ,73°40' ~ 135°5'E) , 
主要 指 淮河 一 秦岭 及 其 西 线 的 延伸 及 其 以 北 地 区 ， 
土地 面积 约 5.78 x 10° km 。 区 内 地 貌 形态 复杂 多 
样 ,包括 山地 、 高 原 、 倪 地 、 丘 陵 、 平 原 等 。 地 势 西高 
东 低 ,形成 层 层 降 低 的 阶梯 状 斜 面 。 研 究 区 主要 位 
于 温带 气候 区 域 , 自 东 向 西 因 水 分 递减 而 分 为 湿润 、 
亚 湿润 、 亚 干旱 以 及 干旱 气候 区 。 区 内 植被 类 型 多 
样 ,形成 了 多 样 而 独特 的 生物 群落 。 
1.2 研究 数据 

遥感 数据 来 自 最 新 发 布 的 GIMMS NDVI 3g v 
1.0 数据 集 ,时 间 跨 度 自 1981 年 7 月 至 2015 年 12 
月 。 数 据 下 载 网址 为 https://ecocast. arc. nasa. gov/ 
data/pub/gimms/3g. v1/。 基 于 Matlab 软件 对 数据 
进行 预 处 理 ,将 全 球 原始 的 ncd 格式 数据 转 成 tf 格 
式 的 中 国 北方 NDVI 数据 集 。 采 用 最 大 值 合成 法 
(MVC) °° ,获得 月 NDVI 值 ,对 月 NDVI 求 平均 得 
到 生长 季 NDVI。 本 研究 将 NDVI <0. 1 的 地 区 认为 
是 “无 植被 区 ” ,不 做 分 析 ” 1。 考 上 处 中 国 北 方 不 同 
地 区 植物 生长 季 始 期 与 末期 不 同 , 且 中 国 北方 冬季 
大 部 分 植物 停止 生长 或 被 冰雪 覆盖 ,所 以 ,确定 5 一 
9 月 为 植物 生长 季 '” 。 
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图 1 研究 区 植被 类 型 


Fig.1 Main vegetation forms in the study area 


气象 数据 来 自 中 国 气象 科学 数据 共享 服务 网 
(http ://cde. ema. gov. en) ,选取 中 国 北方 15 省 共 
408 个 气象 站 点 1982 一 2015 年 逐日 气象 资料 ,计算 
极端 气温 指数 。 

植被 数据 来 自 中 国 科学 院 资源 环境 科学 数据 中 
iÒ ( http :// www. resde. cn) 。 为 便于 分 析 , 将 相似 植 
被 类 型 合并 为 一 类 ,最 终 将 中 国 北方 植被 类 型 分 为 
针 叶 林 、 针 阀 混 交 林 、 落 叶 阔 叶 林 、 灌 从 、 充 漠 植 被 、 
草地 、 高 山 植被 和 栽培 植被 等 8 类 (图 1)。 

1.3 研究 方法 

1.3.1 极端 气温 指数 ”本 研究 选取 气候 变化 检测 
和 指数 专家 组 (expert team on climate change detec- 
tion and indices ) 推荐 的 极端 气温 指数 ( 表 1)。 本 研 
究 将 极端 气温 指数 分 为 冷 极 值 指数 、 暧 极 值 指 数 以 
及 其 他 指数 3 类 。 采 用 基于 R 语言 编辑 的 RClim- 
Dex 程序 完成 指数 计算 。 

1.3.2 极点 对 称 模 态 分 解 方法 极点 对 称 模 态 分 
解 方 法 (extreme-point symmetric mode decomposition， 
ESMD ) 是 对 经 验 模 态 分 解 方 法 (empirical mode de- 
composition ,EMD ) 的 改进 ,该 方法 擅长 寻找 数据 的 
变化 趋势 ,能够 从 观测 序列 中 分 离 出 年 际 变化 趋势 
和 气候 变化 总 趋势 ,是 目前 提取 信号 变化 趋势 的 最 
新 方法 之 一 " 。ESMD 方法 模 态 分 解 主要 将 原始 
HIREA X AMJ X = YM, +R, BVA HAY ESMD 
将 时 间 序 列 X 分 解 成 一 系列 经 验 模 态 M 和 一 个 趋 
AM Ro 

本 研究 除了 使 用 线性 趋势 方法 分 析 1982 一 
2015 年 中 国 北 方 气温 极 值 和 NDVI 的 变化 趋势 外 ， 
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P 国 北方 生长 季 NDVI 变化 及 其 对 气温 极 值 的 响应 


表 1 极端 气温 指数 介绍 


Tab.1 Introduction of the 18 extreme temperature indices 


指标 类 型 ID 指数 名 称 定义 单位 
冷 极 值 FDO 霜冻 日 数 年 内 日 最 低 气 温 <0 °C 的 日 数 d 
IDO 冰冻 日 数 年 内 日 最 高 气温 <0 °C AY H Be d 
CSDI 冷 持续 日 数 年 内 连续 6 d 最 低 气温 小 于 第 10 个 百 分 位 数值 的 日 数 d 
TNmean 日 最 低 气温 平均 值 日 最 低 气 温 的 平均 值 TC 
TN10p 冷 夜 日 数 年 内 日 最 低 气 温 小 于 第 10 个 百 分 位 数值 的 日 数 d 
TX10p 冷 层 日 数 年 内 日 最 高 气温 小 于 第 10 个 百 分 位 数值 的 日 数 d 
TNn 日 最 低 气 温 极 低 值 每 月 日 最 低 气 温 的 最 小 值 C 
TXn 日 最 高 气温 极 低 值 每 月 日 最 高 气温 的 最 小 值 TC 
暧 极 值 SU25 夏季 日 数 年 内 日 最 高 气温 >25 °C 的 日 数 d 
TR20 热 夜 日 数 年 内 日 最 低 气 温 >20 °C 的 日 数 d 
WSDI 暧 持续 日 数 年 内 连续 6 d 最 高 气温 大 于 第 90 个 百 分 位 数值 的 日 数 d 
TXmean 日 最 高 气温 平均 值 日 最 高 气温 的 平均 值 C 
TN90p 暖 夜 日 数 年 内 日 最 低 气温 大 于 第 90 个 百 分 位 数值 的 日 数 d 
TX90p 暧 层 日 数 年 内 日 最 高 气温 大 于 第 90 个 百 分 位 数值 的 日 数 d 
TNx 日 最 低 气温 极 高 值 每 月 日 最 低 气 温 的 最 大 值 TC 
TXx 日 最 高 气温 极 高 值 每 月 日 最 高 气温 的 最 大 值 GS 
其 他 DTR 气温 日 较 差 年 内 日 最 高 气温 和 日 最 低 气温 的 差 值 C 
GSL 生长 季 日 数 连续 6 d 平均 气温 >5 及 连续 6 d 平均 气温 <5 C 的 时 间 蜂 度 d 


还 利用 ESMD 方法 对 气温 极 值 和 NDVI 序列 提取 分 
析 过 去 34 a 的 长 期 变化 趋势 。 

1.3.3 其 他 方法 “Mann-Kendall 法 作为 一 种 非 参 
数 检 验方 法 ,直接 采用 观测 数据 的 秩序 , 而 不 必 事 先 
假定 数据 的 分 布 特征 ,能 较 好 地 处 理 缺 失 值 和 异常 
值 ,适用 于 气象 水文 数据 时 间 序 列 变 化 的 趋势 分 
BO ,本 研究 采用 Mann-Kendall 趋势 检验 方法 
对 气候 极 值 的 变化 趋势 进行 显著 性 检验 。 此 外 ,一 
元 线性 回归 趋势 分 析 法 被 用 于 模拟 每 个 栅 格 的 变化 
趋势 ,反映 NDVI 变化 趋势 的 空间 特征 ~ 
son 相关 系数 被 用 于 量化 NDVI 与 不 同 极端 气温 指 
数 之 间 的 相关 性 。 


> Pear- 


2 结果 与 分 析 
2.1 植被 覆盖 变 
图 2 显示 了 近 34 a 中 国 北方 及 不 同 植被 类 型 


生长 季 NDVI 的 时 间 变 化 特征 。 结 果 表 明 , 中 国 北 
方 生 长 季 NDVI 以 0. 002 . (10a) 一 的 速率 增加 
(P <0.05) ,植被 覆盖 状况 改善 。ESMD 方法 显示 ， 
1982—1992 年 NDVI 逐渐 升 高 ,直至 2005 年 略 有 下 
降 ,随后 又 逐渐 升 高,NDVI 波动 增加 。 

不 同 植被 类 型 的 变化 特征 各 不 相同 。 线 性 趋势 
方法 显示 针 叶 林 灌 丛 .荒漠 植被 .草地 以 及 栽培 植 


被 的 NDVI 呈 增 加 趋势 , 增 速 分 别 为 0. 003 - 
(10a) -0.004. (10a)™!'(P <0.01).0.002. 
(10a)! (P <0. 05) .0.003. (10a) ~' (P <0. 05) VI 
W. 0.008 + (10a) ~! (P <0. 001) ,栽培 植被 增 速 最 
快 。 针 阔 混 交 林 .落叶 阔 叶 林 和 高 山 植 被 呈 减 少 趋 
势 , 减 小 速率 分 别 为 -0.003 . (10a) ~', - 0.001 - 
(10a) ~ FI -0.002 . (10a) ,变化 不 显著 。ESMD 
方法 显示 , 针 叶 林 和 落叶 阔 叶 林 均 表 现 出 先 增加 ， 
1990s 中 期 开始 减少 ,2005 年 后 趋 于 平缓 的 变化 趋 
势 ; 针 阔 混 交 林 HEM .荒漠 植被 草地 和 栽培 植被 的 
变化 趋势 与 整个 中 国 北方 NDVI 的 变化 趋势 较为 相 
似 , 均 表现 出 先 增加 ,1990s 中 期 开始 减少 ,2000s 中 
期 开始 增加 的 趋势 ;高 山 植被 表现 出 先 减少 后 增加 ， 
2000 年 后 呈 缓 慢 的 下 降 趋势 。 

图 3 显示 了 中 国 北方 生长 季 NDVINDVI 变 化 
趋势 以 及 NDVI 显著 变化 的 空间 分 布 情 况 。 
NDVI <0.1 的 地 区 主要 集中 在 西北 地 区 ,包括 新 疆 
中 部 .青海 北部 .甘肃 西北 部 和 内 蒙古 西部 。 中 国 北 
F NDVI 值 从 东南 向 西北 逐渐 降低 , NDVI 高 值 区 主 
要 集中 在 平原 .山地 河谷 和 绿洲 。 黄 土 高 原 地 区 生 
KÆ NDVI 增 速 最 快 ,植被 退化 速度 最 快 的 区 域 主 
要 集中 在 大 小 兴安 岭 以 及 长 白山 一 带 。NDVI 增加 
的 区 域 超过 全 区 62.8% ,显著 增加 的 区 域 超过 
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图 2 1982—2015 年 中 国 北方 不 同 植被 类 型 生长 季 NDVI 的 年 际 变化 及 趋势 


Fig.2 Interannual variations and trends of the NDVI of different vegetation forms in North China in growing season from 1982 to 2015 
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南山 区 .黄土 高 原 .河套 平原 .吕梁 山 和 太行 山大 别 
山 以 及 辽西 丘陵 地 区 。 减 少 的 面积 约 占 整个 地 区 的 
37.2% ,显著 减少 的 面积 约 12% , 主要 分 布 在 大 兴 
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北方 生长 季 NDVI 空间 变化 
Fig.3 Spatial variation of the NDVI in North China in growing season from 1982 to 2015 


BUR ,小兴安岭 和 长 白山 区 。 


2.2 极端 气温 变 


化 


2.2.1 冷 极 值 图 4 显示 了 气温 冷 极 值 指 数 的 时 
间 变 化 趋势 , 表 2 为 Mann-Kendall 趋 势 检验 结果 
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图 4 极端 气温 指数 的 年 际 变化 及 趋势 
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Interannual variations and trends of the extreme temperature indices in North China from 1982 to 2015 
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R2 冷 极 值 气温 指数 趋势 检验 
Tab.2 Test results of the trends of cold extreme 


temperature indices 


指数 趋势 检验 单位 
FDO -4.534 =" d- (10a) 7! 
CSDI -0.459 d- (10a) 7! 
TN10p -2.585 *™ d» (10a) ~! 
TNn 0.154 T+ (10a) 7! 
IDO -1.006 d- (10a) 7! 
TNmean 0.427 *™* c+ (10a) 7! 
TX10p 一 1.783 ** d+ (10a) ~! 
TXn -0.088 T+ (10a) 7! 


注 ; * 表示 相关 性 在 P<0.05 水 平 显著 ; * * 表示 相关 性 在 P<0.01 
水 平 显著 ; * * * 表示 相关 性 在 P<0.001 水 平 显著 。 下 同 。 


结果 显示 ,多 数 冷 极 值 指 数 呈 下 降 趋势 。 近 34 a 
X, FDO, TN1IOp 和 TX10p 分 别 以 4. 534 d > 
(10a) ~', - 2. 585 d + (10a) 7' Ñ - 1. 783 d> 
(10a) 的 速率 显著 下 降 , TNmean LI 0. 427 © >» 
(10a) 的 速率 逐渐 升温 (P <0. 001), ESMD 方法 
显示 冷 极 值 指 数 的 变化 趋势 各 异 : FDO, IDO 和 
TX10p 变化 趋势 相似 , 均 表 现 出 先 迅 速 下 降 、2005 
年 后 平稳 变化 的 趋势 ;CSDI 和 TN10p 则 显示 出 先 
下 降 、1995 年 左右 达到 最 低 值 ,后 波动 上 升 的 变化 
趋势 。 多 数 冷 极 值 指 数 呈 下 降 趋 势 , 说 明 随 着 气候 
变 暧 ,日 最 低 气 温 逐 渐 升 高 ,表征 低温 事件 的 极端 气 
温 冷 极 值 指数 多 数 呈 下 降 趋势 。 

2.2.2 RAIE 图 4 显示 了 气温 暧 极 值 指数 的 时 
间 变 化 趋势 , 表 3 为 Mann-Kendall 趋势 检验 结果 。 
所 有 上 暧 极 值 指数 均 呈 显著 上 升 趋势 , WSDI、TXx、 
SU25 .TR20 .TXmean , TN90p ,TX90p 和 TNx 分 别 以 
1.251 d- (10a) ~! 0.449 C . (10a) 7! 4.704 d- 
(10a) ~' 2.758 d+ (10a) 7',0.34 C + (10a) 7', 
4.32 d+ (10a) ~! 3.183 d+ (10a) “'#l0. 456 C >» 


表 3 暖 极 值 气温 指数 趋势 检验 
Tab.3 Test results of the trends of warm extreme 


temperature indices 


指数 趋势 检验 单位 
SU25 4.704 ** d+ (10a) ~! 
WSDI 1.251 * d+ (10a) 7! 
TN90p 4.32 ** d+ (10a) ~! 
TNx 0.456 ** + (10a) ~! 
TR20 2.758 “* d+ (10a) ~! 
TXmean 0.34 ** T+ (10a) 7! 
TX90p 3. 183 ** d+ (10a) 7! 
TXx 0.449 * C+ (10a) ~! 


(10a) ”的 速率 呈 上 升 趋势 。ESMD 估计 方法 显示 ， 
暖 极 值 指数 大 致 表现 了 自 1982 年 开始 波动 上 升 并 
在 2000s 初 略 有 下 降 的 变化 趋势 。 日 最 高 气温 逐渐 
升 高 ,表征 高 温 事 件 的 极端 气温 暧 极 值 指数 呈 上 升 
趋势 ,对 全 球 气候 变 暧 呈正 响应 。 

2.2.3 ”其 他 指数 ”图 4 显示 了 气温 日 较 差 和 生长 
季 日 数 的 时 间 变 化 趋势 , 表 4 为 Mann-Kendall 趋势 
检验 结果 。DTR 以 -0.087 € + (10a) 一 的 速率 呈 
减 小 趋势 (P <0.05)。GSL 以 4.401 d+ (10a) HY 
速率 呈 增 加 趋势 (P <0. 001), ESMD 估计 方法 显 
示 ,DTR 先 减 小 后 增 大 并 于 1997 年 左右 达到 最 大 
值 , 之 后 逐渐 减 小 ;CSL Á 1982 年 起 逐渐 增 大 ,2005 
年 左右 达到 最 大 值 后 略 有 减 小 。 


表 4 其 他 气温 指数 趋势 检验 
Tab.4 Test results of the trends of other 


temperature indices 


指数 趋势 检验 单位 
DTR -0.087 * + (10a) >! 
GSL 4.401" d+ (10a) -1 


2.3 NDVI 对 极端 气温 指数 的 响应 

图 5 表示 中 国 北方 及 NDVI 显著 变化 区 域 对 极 
端 气 温 指 数 的 响应 。NDVI 显著 增加 区 域 主要 包括 
天 山 塔里木 盆地 北 缘 ,祁连山 区 .陇南 山地 河套 平 
原 .黄土 高 原 ,太行山 和 吕梁 山区 .辽西 丘陵 .大别山 
区 ;显著 减少 区 域 主要 包括 大 兴安 岭 小兴安岭 K 
白山 。 
2.3.1 对 冷 极 值 的 响应 气温 冷 极 值 指数 中 ， 
FDO IDO .CSDI TN10p 及 TX10p 与 中 国 北方 NDVI 
呈 负 相关 ,其 中 FD0 .TN10p #1 TX10p 与 NDVI 的 相 
关系 数 通过 95% 显著 性 检验 (图 5) FDO 和 TN10p 
代表 夜间 的 低温 状态 ,TX10p 代表 白天 的 低温 状态 ， 
气温 过 低 易 导致 积温 不 足 ,严重 时 引发 低温 冻 灾 , 影 
响 植物 生长 。 负 相关 说 明 随 着 气候 变 暧 ,霜冻 风险 
降低 ,夜间 变 暧 为 植物 生长 创造 了 更 适宜 的 条 件 ; 冷 
昼 日 数 减少 ,日 间 气温 升 高 更 有 利于 植物 光合 作用 ， 
保障 植物 生长 。TNmean ,TNn 和 TXn 与 NDVI 呈正 
相关 ,其 中 TNmean 与 NDVI 的 相关 系数 通过 99% 
显著 性 检验 。TNmean 代表 白天 的 低温 状态 ,日 最 
低温 升 高 ,使 植被 受 日 间 低温 伤害 的 风险 降低 ,植被 
光合 作用 加 强 ,保障 并 促进 了 植物 生长 发 育 。 

各 地 区 对 冷 极 值 指数 的 响应 程度 不 同 。NDVI 
显著 增加 区 域 与 多 数 冷 极 值 指数 以 负 相 关 为 主 ， 
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注 : 显 著 相关 指 在 95% 水 平 上 通过 显著 性 检验 , 极 显 著 相关 指 在 99% 水 平 上 通过 显著 性 检验 。 


图 5 中 国 北方 及 不 同 区 域 极 端 气温 指数 与 NDVI 的 相关 性 统计 


Fig.5 The correlation coefficients between the NDVI and the extreme/temperature indices in North China and different regions 


FDO 和 TNmean 是 NDVI 增加 区 响应 强烈 的 指数 。 
塔里木 岔 地 北 缘 .河套 平原 太行山 和 上 吕梁 山区 等 对 
冷 极 值 指数 响应 强烈 。NDVI 显著 减少 区 域 与 多 数 
冷 极 值 指数 呈正 相关 ,但 响应 情况 各 异 : 大 兴安 岭 区 
对 TNn 响应 强烈 ,小兴安岭 区 对 TN10p 响应 强烈 ， 
长 白山 区 则 对 TX10p 响应 强烈 。 
2.3.2 对 暖 极 值 的 响应 8 个 气温 暖 极 值 指数 和 
中 国 北方 NDVI 均 呈 正 相 关 , TR20, TXmean 及 
TX90p 与 NDVI 的 相关 系数 通过 95% 显著 性 检验 ， 
TN90p 与 NDVI 的 相关 系数 通过 99% 显著 性 检验 
(图 5)。 说 明 随 着 气温 升 高 ,适宜 植物 生长 的 暧 区 
范围 扩展 ,加 上 夜间 温度 升 高 .温暖 夜 数 增多 有 效 减 
少 了 植被 受到 极端 低温 伤害 的 风险 ,有 利于 植物 的 
生长 发 育 。 

从 整体 来 看 , NDVI 显著 增加 区 域 对 暖 极 值 指 
数 以 正 响应 为 主 , 显 著 减 少 区 域 则 多 呈 负 响应 ， 
NDVI 显 著 增 加 区 域 对 TN90p 以 显著 正 响应 为 主 ， 
NDVI 显著 减少 的 区 域 对 SU25 呈 显 著 负 响应 (P < 
0.05), SU25 可 表征 为 夏季 日 间 的 高 温 状态 ,NDVI 
显著 减少 区 域 主要 分 布 在 东北 林 区 ,生长 季 长 度 在 
120 ~ 160 d" ,夏季 高 温 日 数 多 、 高 温 持 续 时 间 长 
易 破 坏 森 林 植 物 生长 环境 的 平衡 性 ,近年 来 
SU25 持续 上 升 ,可 能 造成 森林 植物 生长 受 损 ,生物 
量 减少 ,植被 覆盖 退化 。 天 山 .塔里木 盆地 北 缘 、 祁 
连 山 区 ,太行山 和 吕梁 山区 .小 兴安 岭 等 区 域 对 暖 极 
值 指数 响应 强烈 。 
2.3.3 对 其 他 指数 的 响应 ”图 5 显示 中 国 北方 


NDVI 与 DTR 呈 负 相关 与 GSL 呈正 相关 , CSL 与 
NDVI 的 相关 系数 通过 95% 显著 性 检验 。 说 明 随 着 
气温 升 高 ,生长 季 长 度 延 长 ,植物 生长 时 间 更 加 充 
分 ,更 有 益 于 植被 的 生长 发 育 。 

DTR 指 日 间 气温 和 夜间 气温 的 差 值 , 可 表征 为 
促进 植物 生长 发 育 的 有 效 温度 。DTR 在 多 数 区 
域 主要 以 负 相 关 为 主 ,天 山 、 辽 西 丘陵 、 大 兴安 岭 尤 
为 显著 ,主要 是 由 于 DTR 呈 减 小 趋势 ,气温 日 较 差 
缩小 不 利于 有 机 质 的 积累 , 随 着 夜间 温度 升 高 ,植被 
呼吸 消耗 增加 ,有 机 质 积累 量 减少 ,植物 生长 受到 限 
制 。GSL 在 大 部 分 区 域 呈正 相关 , 且 相 关 性 强 的 区 
域 主要 集中 于 NDVI 增加 区 ,进一步 说 明生 长 季 延 
长 对 植物 生长 有 积极 作用 。 


3 ”结论 与 讨论 


3.1 结论 

(1) NDVI 的 时 间 变 化 显示 ,1982 一 2015 年 中 
国 北方 生长 季 NDVI 以 0.002. (10a) -的 速率 增加 
(P <0. 05 ) ;ESMD 方法 表明 ,NDVI 波动 增加 。 针 
叶 林 、 灌 从 、 荡 漠 植被 .草地 以 及 栽培 植被 显 车 增加， 
其 中 栽培 植被 增 速 最 快 ; 针 冰 混交 林 落叶 阔 叶 林 和 
高 山 植被 呈 不 显著 减少 趋势 。 

(2) NDVI 的 空间 变化 显示 ,黄土 高 原 地 区 生 
长 季 NDVI 增 速 最 快 ,NDVI 退化 速度 最 快 的 区 域 主 
要 集中 在 大 、 小 兴安 岭 及 长 白山 。NDVI 显著 增加 
的 区 域 超过 全 区 的 33% ,主要 分 布 在 天 山 、 塔 里 木 
盆地 北部 、 祁 连 山 、 陀 南山 区 、 黄 土 高 原 、 河 套 平原 、 
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吕梁 山 和 太行 山 、 大 别 山 以 及 辽西 丘陵 地 区 。 显 著 
减少 的 面积 约 12% ,主要 分 布 在 大 兴安 岭 . 小 兴安 
岭 和 长 白山 区 。 

(3) 气温 极 值 变 化 显示 , 除 日 最 低 气温 平均 值 
(TNmean) 和 日 最 低 气 温 的 极 低 值 (TNn) 呈 上 升 趋 
势 外 ,其 余 冷 极 值 指数 都 表现 为 下 降 趋势 ; 所 有 了 暧 极 
值 指数 均 呈 上 升 趋势 ;其 他 指数 中 ,DTR StF BERS 
势 ,GSL 呈 上 升 趋势 。 

(4) 从 NDVI 与 极端 气温 指数 的 相关 性 看 , 冷 
极 值 指数 中 FD0 IDO .CSDI、TN10p 及 TX10p 与 ND- 
VI 2 fi 4H% ,TNmean TNn 和 TXn 45 NDVI 呈正 相 
3, FDO TN10p .TX10p 及 TNmean 的 相关 系数 通过 
显著 性 检验 ;所 有 了 暧 极 值 指数 与 NDVI 均 呈 正 相 关 ， 
TR20 .TXmean .TX90p 及 TN90p 的 相关 系数 通过 显 
著 性 检验 ;NDVI 与 DTR 呈 负 相关 ,与 GSL 呈正 相 
关 ,GSL 与 NDVI 的 相关 系数 通过 显著 性 检验 。 

(5) RU 塔里木 盆地 北 缘 .祁连山 区 、 河 套 平 
原 黄土 高 原 、 太 行 山 和 吕梁 山区 等 NDVI 显著 增加 
区 域 对 极端 气温 指数 的 响应 强烈 。NDVI 显著 增加 
区 对 FDO \TNmean ,TN90p 和 GSL 等 指数 响应 较 强 。 
NDVI 显著 减少 的 大 、 小 兴安 岭 及 长 白山 区 对 极端 
气温 指数 的 响应 各 异 , 主要 与 SU25 呈 显 著 负 相关 
(P <0.05), 

3.2 讨论 

气温 极 值 对 植被 的 影响 是 一 个 复杂 的 过 程 。 本 
研究 分 析 了 NDVI 对 全 部 极端 气温 指数 的 响应 情 
况 , 旨 在 确定 对 植物 生长 影响 较 强 的 极端 气温 指数 ， 
以 便 未 来 分 析 区 域 植物 生长 对 极端 气候 的 响应 时 可 
以 选择 更 有 指示 意义 的 指数 。 中 国 北 方 NDVI 与 
FDO, TNmean, TN10p, TX10p, TR20, TXmean , 
TX90p .TN9Op 及 GSL 存在 强 相关 性 ,多 数 NDVI W 
著 增 加 的 区 域 对 FDO TNmean ,TN90p 和 GSL 响应 
强烈 。Yao 等 中 研究 也 指出 ,新 疆 地 区 NDVI 的 变 
化 与 TNmean 呈正 相关 ,与 TN90p 呈 负 相关 ; 赵 安 周 
等 (中 研究 发 现 , 陕 甘 宁 地 区 NDVI 与 TXmean、 
TNmean 存在 显著 相关 性 。 未 来 研究 中 国 北方 植被 
对 极端 气温 的 响应 时 ,可 以 优先 考虑 以 上 与 NDVI 
存在 强 相关 的 极端 气温 指数 。 但 也 有 研究 结果 显 
示 , 郡 阳 河 流域 TXn TNn, TNx 和 TXx 对 NDVI 的 影 
响 较 强 ”” ;福建 省 TN90p、TX90p 对 植物 生长 的 影 
MER 。 可 见 , 植 被 覆盖 对 极端 气温 的 响应 是 复 
杂 而 多 样 的 。 目 前 对 这 两 者 之 间 的 响应 缺乏 大 量 系 
统 的 分 析 和 研究 。 本 研究 揭示 了 中 国 北方 植被 对 极 


端 气温 的 响应 情况 ,对 植被 覆盖 与 极端 气候 响应 研 
究 进行 了 补充 。 但 并 未 对 植被 覆盖 对 极端 气温 指数 
的 响应 在 空间 上 的 表现 情况 做 深入 研究 ,未 来 需要 
进一步 分 析 植 被 变化 对 极端 气温 指数 响应 的 空间 表 
现状 况 。 另 外 ,本 研究 仅 分 析 了 极端 气温 对 NDVI 
的 影响 ,未 将 极端 降水 变化 对 植物 生长 的 影响 纳入 
研究 范围 。 未 来 在 研究 极端 气候 事件 与 植被 覆盖 之 
间 的 相关 关系 时 ,应 考虑 气温 和 降水 之 间 的 协同 效 
应 ,探究 植物 生长 与 各 因子 之 间 的 机 理 模 型 ,这样 就 
能 定量 分 析 各 影响 因子 对 植物 生长 的 贡献 率 ,更 准 
确 地 识别 极端 气候 条 件 下 中 国 北方 生态 系统 的 脆弱 
性 和 复杂 性 。 
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Spatiotemporal Change of NDVI and Its Response to Extreme 
Temperature Indices in North China from 1982 to 2015 


HE Hang, ZHANG Bo, HOU Qi, LI Shuai, MA Bin, MA Shang-gian 
(College of Geographical and Environmental Sciences , Northwest Normal University , Lanzhou 730070 , Gansu , China ) 


Abstract: Climate warming is conducive to enhancing vegetation activities. Here ,the interannual and spatial vari- 
ations of vegetation cover in north China in growing season were analyzed based on the satellite-derived normalized 
difference vegetation index ( NDVI) , and its responses to the change of extreme temperature indices were studied by 
using GIMMS NDVI 3g V1.0 datasets, daily temperature and precipitation data. Across the whole study area, the 
trends calculated by linear regression showed that the NDVI value in growing season increased at a rate of 0.002 + 
(10a) ~' from 1982 to 2015. The results of extreme-point symmetric mode decomposition showed that the NDVI in- 
creased gradually until 1992 , decreased slightly until 2005 , and then increased gradually. The NDVI values of conif- 
erous forest , shrubbery , desert vegetation , grassland and cultivated vegetation were all in an increase trend, and those 
of mixed forest, broadleaved deciduous forest and alpine vegetation were in a decrease trend. Spatially, the NDVI 
was in a decrease trend from the southeast to the northwest, and the area of the regions where the vegetation was sig- 
nificantly improved accounted for 33% of north China. The regions where the NDVI increased significantly were 
mainly distributed in the Tianshan Mountains and north Tarim Basin in north Xinjiang, Qilian Mountains , mountain- 
ous area in south Gansu Province, Loess Plateau , Hetao Plain, Lvliang Mountain , Taihang Mountain, and hilly region 
in west Liaoning Province. The area of the regions where the NDVI decreased significantly were mainly distributed in 
the Great Khingan Range, Lesser Khingan Mountains and Changbai Mountain. Among the 18 extreme temperature 
indices ,except the mean daily minimum air temperature and the lowest minimum air temperature were in an in- 
crease trend, all others of cold extreme temperature indices were in a decreased trend; the warm extreme tempera- 
ture indices were all in an increase trend. The NDVI was negatively correlated with FDO ,TN10p and TX10p (P < 
0.05) ,but positively correlated with TNmean (P <0.01).The NDVI was positively correlated with all warm ex- 
treme temperature indices, and was significantly correlated with TR20,TXmean, TX90p and TN90p (P <0.05). 
There was also a significant positive correlation between NDVI and GSL (P <0.05). 

Key words: growing season; GIMMS NDVI 3g V1.0; spatiotemporal variation; extreme temperature index; re- 


sponse; North China 


